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I. INTRODUCTION 


A. MID-ATLANTIC BIGHT OBJECTIVES 

Better knowledge of operating environments may allow any Navy to improve the 
use of its USW systems in different tactical scenarios. Even though undersea sound and its 
practical uses have been continuously explored since World War II, testing of other 
techniques are necessary in order to obtain important sea environment information which 
may help in the development of accurate Naval USW systems. To face this challenge, the 
Mid-Atlantic Bight (MAB) experiment was conducted. The MAB experiment was an 
intensive field study of the shelfbreak south of Long Island dunng the strongly stratified 
summer season whose basic objectives were: 

1. To provide detailed time series of the oceanographic conditions along typical 
major acoustic paths by using moored acoustical arrays. 

2. To provide a diversity of frequencies and source locations from broadband SUS 
sources. 

3. To determine the temperature structure of the ocean environment by using 
acoustic tomography. 

This study may provide the Navy an alternative and simple means of determining 
characteristics in unknown or poorly known operating environments. (Pickart et al., 1996) 
B. STUDY OBJECTIVES 

The objectives for reaching the fundamental goal of this study are as follows: 

1. To determine the bottom sound speed by measuring the coherence of two 
vertically displaced sensors in the shallow water channel and comparing this 
function with one obtained from a simple model equation over a range of 
frequencies. 

2. To determine the bottom density by comparing mode travel times of computed 
spectrograms for a phone from SUS explosive data and those obtained by 


employing the concept of an effective depth for an ideal waveguide. 


] 


3. To determine the bottom attenuation by finding mode attenuation as a function 


of range extracted from spectrograms and relating it to loss per bottom bounce. 


Il. EXPERIMENTAL/ENVIRONMENTAL DESCRIPTION 


The summer Mid-Atlantic Bight experiment consisted of a cruise from 19 July to 9 
August 1996 aboard the R/V Endeavor from the University of Rhode Island. The cruise 
included three major components: acoustic tomography, SEASOAR operations, and 
hydrography/tracers. An AXBT flight was also conducted during this time period, as well 
as a SUS drop fly-over. The acoustic portion of the experiment consisted of three 
components: 

1. The moored acoustic transmission array, or tomography array 

.2. The moored physical oceanography array, which provided detailed time series 
of the oceanographic conditions along our major acoustic path 

3. The broadband, explosive (MK61 SUS shot) source component. 

Multiple lines of SUS drops were deployed in the area by P-3 aircraft and monitored by two 
vertical line arrays (VLAs) moored on the shelf in ~85m water. (Pickart et al., 1996) 

To achieve the objectives of this thesis, eleven broadband SUS shot transmissions 
recorded at the NW VLA were selected. All of these were dropped along a single track 
from a point at latitude 40°16.2’N, longitude 71°10.2’W to a point at latitude 40°16.8’°N, 
longitude 70°57.0’W. The sampling rate was 3906.25 Hz, the duration of one record was 
8.3886 sec and the number of samples/record was 32768. The arrangement and the location 
of the acoustical elements used for this study are shown in Figure 1. This region of the 
experimental site has a nominal water depth of ~ 90 m with little variation (<10m) from the 
VLA to the positions of the selected SUS drops. The receiving VLA had 16 hydrophones 
of which only the upper 8 successfully recorded the SUS transmissions with an upper depth 
of 30.5 m and inter-element spacing of 3.5m. A typical sound speed profile measured near 


the trial site and used in this analysis is shown in Figure 2. 
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Figure 1. Arrangement and Location of Acoustic Elements Used at the Experimental Site. 
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Figure 2. Sound Speed Profile Measured Near the Trial Site and Used in This Study. 





Ill. ANALYSIS METHODS 


A. INVERSION OF BOTTOM SOUND SPEED 

In determining the bottom sound speed a method was performed to determine the 
critical angle of the seabed from the vertical directionality of the ambient noise in the water 
column (Buckingham and Jones, 1987). The properties of the ambient noise that we are 
specifically concerned with are the vertical directionality and the spatial coherence of the 
fluctuations of two sensors vertically displaced in the noise field. These two quantities are 


related through a Fourier transform relationship, which is given by 


(a) => [ F()explcecos(6)] sin()26 (1) 


where 9 is the polar angle measured down from the zenith and F(6) is a dimensionless 
directional density function which represents the distribution of the noise power in the 


vertical. F(8) can be normalized according to the condition 
I 
5 | F(@)sin(6)a0 = 1. (2) 


This indicates that in the limit as Q —> 0 the value of the coherence function goes to unity. 
Buckingham and Jones provided a simple example for isotropic noise in which F(@) 


is independent of 8. According to Eqs. (1) and (2) the coherence function in this case is 


shown to be 
sin(©2) 
lr... (€2) = ———,, 3 
ee | ) Q ( ) 
where {2 1s a dimensionless angular frequency defined as 
ol 
ae (4) 
e 


Ww 


® is the angular frequency, c, is the sound speed in the water column, and / is the 


Ww 


separation between the two sensor positions. This example is shown in Figure 3. 











Figure 3. Example of Ideal Vertical Coherence Function Versus Frequency. 


This function predicts that the coherence of the noise at two sensors is real due to 
symmetry of the noise field about the horizontal. For an asymmetrical noise distribution in 
the vertical, the coherence function is complex with a large imaginary component. A 
symmetrical noise distribution is also predicted by Buckingham and Jones’ theoretical 
model of shallow water noise, where the coherence function for shallow water noise is an 
oscillatory function that decays with increasing (2. At a given frequency, the directional 


density function of the noise predicted by the model is closely approximated by 


a, for 0, |s0s{ +0, 
F,(9)= 2 2 ; (5) 


0, otherwise 


where a, is independent of 8. Then, the noise propagates in directions close to the 


horizontal, with no acoustic energy traveling with grazing angles greater than @ (the critical 
grazing angle). The contribution of local, wind-driven sources to the overall noise field is 
approximately isotropic, that is, the directional density function of the continuous 
component of the noise can be expressed as 


F,(0) =a,, for0<O<xz, (6) 


where a, is independent of 8. By superposing Eqs. (5) and (6) we have a representative 
directional density function for the total noise field in shallow water: 


F,,(0)= F,0)+ ©) 


7U 7U 
|teq+a) fof =-0, }<0<{ = 40, | (7) 


a otherwise. 


Figure 4 shows a sketch of the vertical distribution expected in a shallow water 
environment. Through the normalization condition in Eq. (2), the critical grazing angle can 
be expressed in terms of the other two parameters as 


laa 


sin(0, ) = . (8) 


a, 





and the model coherence function of the noise represented by the distribution in Eq. (7) is 
found from Eq. (1) to be 


r,,(8) = -- sin(Q) + ante) sin] Q sin(0,.)) ; (9) 


Using this model, the spatial coherence of two vertically displaced sensors of the 


NW VLA was computed from extracted ambient noise data. By spatial coherence is meant 





Figure 4. Model of the Vertical Distribution of Ambient Noise in a Pekeris Waveguide. 
9 


the time-averaged crosscorrelation coefficient of the noise observed by hydrophones 


separated vertically in the water column, 1.e., 
P(e) 


Lara = lim (Gare yee): (10) 


xx ay 


where P,(w) represents the cross spectral density as a function of frequency for two 
separated phones and P,(o) and P,(w) represent the spectral density for each one of the 


separated phones. To search the two unknown parameters a,and 9.in Eq. (9), an 
optimization procedure was implemented by finding a least squares fit between the real part 
of the measured coherence function, I,,., and the theoretical coherence function, I’... By 


meas ° S 


defining the ambiguity function 


N 
= 


n=l 


- (11) 





(ol aee.) Re al{T neas($2)} 





values for a, and 9. are obtained when ® is a minimum. 
B. INVERSION OF BOTTOM DENSITY 

To determine the bottom density, the concept of an effective depth for an ideal 
waveguide was employed (Chapman and Ward, 1989). By considering the phase change of 
a plane wave reflected at a fluid/fluid boundary, Weston introduced the concept of the 
effective boundary depth, arguing that one could view the reflection as taking place at an 
imaginary pressure - release boundary located at a specific depth below the true boundary. 
He proposed that the normal modes of the Pekeris model should be similar to those of an 
ideal pressure - release waveguide whose depth was equal to the true water depth plus the 
effective boundary depth. The normal mode wavenumbers for the Pekeris model can then 
be estimated easily without numerical iteration. 

The Pekeris model consists of a uniform water layer of constant depth bounded 
above by a pressure - release surface and below by a uniform fluid half space having 


different density and sound speed. Consider a plane wave of frequency fin a fluid medium 


of sound speed c,, incident upon a plane boundary with another acoustic medium at angle 0 


relative to grazing, as depicted in Figure 5. 


Q B 


Real Boundary 


H | 


= Imaginary Boundary yp 
Figure 5. Geometry Used to Derive The Effective Boundary Depth. 


The ray segments s and s’ are given by 





AH 
s=— (12) 
sin9 
and 
2AH cos’ 8 
a (13) 
sin8 
and the phase of the reflected plane wave at points A and B is given by 
d,=k,s'+¥ (14) 
and 
dO, =2k s+", (15) 


0) 
in which k,, = me Requiring , =, the effective depth of the boundary is defined by 
Seas 


Hf \0)-H =———_.. 
(6) 2k, sin® 


(16) 
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For an isospeed water column with sound speed c, and density p, overlying a 
homogeneous solid medium with compressional sound speed c,, shear speed c,, and 
density »,, Chapman and Ward (1984), have shown (assuming c, < c,,and bottom 
attenuation 1s negligible) that the effective depth is 


2 
C, BO 
H, =H+\1—2) — pa sin@ (17) 
C k 


where H is the water depth and 9. is the grazing critical angle. If c, <<c,, then Eq. (17) is 


reduced to 


a (18) 
kK, sin8. 


The boundary conditions for this ideal waveguide require the pressure to vanish at 


the upper and lower interfaces. Correspondingly, eigenfunctions, w(z), must satisfy 


(Medwin and Clay, 1998) 





w(z)|,_.= 6 and v2) =0. (19) 


Solutions of the eigenfunctions for these boundary conditions are 


NNZ 
z)=sin} —— |, 20 
w(z) | (20) 
and the vertical component of the wavenumber, y, must then satisfy the quantization 
condition 
nt 
=—., 21 
Y= F (21) 


where 7 is the integer and designates the mode number. The honizontal values of the 


wavenumber K, are then given by 


K,=(k-y7)" , (22) 


12 


and the requirement that K, be real gives the mode cut-off. By substituting y, from Eq. 


(21) and k, = se , K, can be written as 
C 


c-2) Ge) a 
oe C, H, ; 


Furthermore, the group velocity of the n™ mode is given by 


eso (24) 


therefore, 


(25) 
Pe: Ea 
of, 
The group velocity u, and the mode travel times are related by 
t= zs ; (26) 
u 


where R is the horizontal distance between source and receiver and t,. is the travel time of 
the n™ mode. Combining Eas. (25) in (26) the mode travel times can be written 


f. Sr ee (27) 


ie 20 ace z 
2 fH, 


From the definition of effective depth, this becomes 





-1/2 


(28) 


Cc 


Ww 


R nisin 8 
AD -— (2nfHisind \ a 
| : 8). 


which is the final expression for the mode travel times . Since we know the real depth H, 
and c,and 9. were determined from the previous analysis, then we can invert for g by 


lie 


measuring the mode travel times over a range of frequencies for different SUS 


transmissions. Assuming p, ~ 1, then g ~ p,. 


&. INVERSION OF BOTTOM ATTENUATION 
Mode attenuation as a function of range can be estimated from several computed 


SUS spectrograms. The mode attenuation per km, a, , satisfies the relation 
ATL, =a, AR, (29) 


where 47Ly, and AR are the mode transmission loss due to attenuation and range between 
SUS drops, respectively. Neglecting volume attenuation in the water column, the decrease 
in mode amplitude is expected to be caused only by bottom attenuation and cylindrical 
spreading (ZL ,,2.4 = 10logR). Since mode n can be thought of as due to rays propagating at 
angle 8, the mode attenuation can be considered the result of multiple bottom reflection 
losses. 

The number of bottom bounces over | km of honzontal propagation can be found 
by considering the angle of elevation or depression @ of the local direction of propagation 
of the waveform. Referring to Figure 6, the angle 8 and the diagonal distance for one 
bounce L, are related by 

sin§ = cial (30) 
L, 
where 7 is the water depth. Defining the total diagonal distance Z corresponding to R = 1 
km of horizontal propagation by 


R 
cos0- 





L= (31) 


then L and L, can be related in order to define the number of bottom bounces per km as 





va £ = Rtane 


32 
Leaps te 


The attenuation for mode 7 can then be expected to be due to N (0, ) bottom bounces per 


km producing 


—_ N(0,)or,,(0,) (33) 


where @_ is the corresponding propagation angle for mode n and a,, (0, )is the bottom loss 


per bounce for a single reflection at this angle. 
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L SS | Bottom 


L 


Figure 6. Geometry Used to Find the Number of Bottom Bounces. 
The reflection coefficient at a fluid-fluid interface for a plane wave traveling at 
angle 0 is defined by 


eo Pro, Sue its IPPs , (34) 
p,¢, sin9 —ip_b, 


where c,, and p, are the sound speed and density, respectively, in the upper medium, c, and 


p,, are the sound speed and density, respectively, in the lower medium, and 


1/2 


D = (2) cos’ 0- ! . (35) 
Cy 


is} 


If there 1s no attenuation, then |R| = 1 for incidence beyond the critical angle. However, if 


the bottom is lossy then 


|R(0) =1-<(6), (36) 
and the bottom loss per bounce can be defined as 
ot, (0,) = -20log(R(6, }) = 20e(, ). (37) 


Combining the equations above, the reduction in the reflection coefficient amplitude can be 
determined by 


1 a, 1 ATL 
(,) : 


~ 20 N(@,)  20ARNO,) w 


Bottom attenuation can be introduced by defining an imaginary component of the 


bottom sound speed, 
©. XR SIS. (39) 


such that the wavenumber in the lower medium becomes 


Ke ee a 1 in.,. (40) 
Cp 
where 
3) 
ques (41) 
Cy 
Since the solution in the lower medium goes like 
Bren = Blane (42) 


then a, is found to produce the exponential decay in range. This then causes transmission 
loss according to 


ATL se, = 8.6860,AR=QAR, (43) 


atten 


where « 1s the bottom attenuation in units dB/m, and can be written as 


a= Bese (44) 
Cy 
In order to estimate the bottom attenuation (dB/m), a simple matching algonthm 


may be employed which searches over values of 6 (the complex bottom sound speed) until 
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the reduction in the reflection coefficient amplitude agrees with the observed value of €(8,) 
defined in Eq. 37. Once values for 6 are obtained (for various frequencies and mode angles 
8,), then Eq. 43 may be employed to produce the final estimate for the bottom attenuation 


x: 


li 





IV. DATA ANALYSIS AND RESULTS 


A. INVERSION OF BOTTOM SOUND SPEED 

As described in Section II. A., the inversion for bottom sound speed is based on the 
evaluation of the vertical coherence of the ambient noise field. For this analysis, the 
ambient noise was extracted from segments of the SUS data sets between actual SUS 
signals (and allowing for reverberation to fade). The spectrum of the ambient noise 
computed showed little variation from low frequency up to 600 Hz indicting no strong 
unexpected noise sources were present. Also, the imaginary component of the coherence 
between the two hydrophones was found to be small as required. This indicates that the 
ambient noise field at the tnal site for frequencies below 600 Hz has a symmetrical 
distribution in the vertical direction and the data is suitable for inversion by the method 
described previously. 

It should be noted that the inversion method described in Section IH. A. was 
originally defined for an isospeed water column overlying a homogenous bottom. As is 
obvious from the local sound speed data (refer to Figure 2), there is significant variation in 
sound speed over depth. However, if the sound speed change is not too large over the 
vertical separation of the phones, the inversion method should still be applicable by simply 
adjusting the local angles to the appropriate values at the bottom interface using Snell’s law. 

Since the data from the upper two phones of the VLA was analyzed, a local sound 
speed of c, ~ 1485 m/s was used in the calculation of the ambiguity function (refer to Eq. 
LOE 

The fit between the theoretical coherence obtained from the model and the real 
component of the observed coherence is shown in Figure 7. By computing the ambiguity 


surface for the optimization of two parameters (a, and 8.) an optimal value was found at 


O.~18° in the water column. Figure 8 shows the ambiguity surface for the optimization 


with two parameters, a, and 6, whose minimum corresponds to a, ~ 0.18 and mean critical 
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angle 0.~+18°+1° for 4 different SUS explosive data. A grazing critical angle of 
approximately 19°+1° near the bottom was determined from Snell’s law. The sound speed 
close to the bottom was 1476 m/s, and therefore the bottom sound speed c,+1563+10 m/s. 
This value was compared to values obtained by a sophisticated genetic algorithm technique 
employed at the University of Rhode Island (Miller and Potty, 1998). They found values in 
the upper few meters of the sediment in agreement with the result reported here. 
B. INVERSION OF BOTTOM DENSITY 

To invert for bottom density, the method of ideal waveguide effective depth was 
introduced. However, this method assumed a constant sound speed in the water column. In 
fact, this was fundamental in the definition of the ideal waveguide eigenfunctions. For that 
reason, this method may not be expected to produce very good results. We will attempt to 
invert for density anyway by assuming the water column is isospeed with nominal sound 
speed c,, ~ 1485 m/s. Consistent with this value, we shall assume the critical angle 0, ~ 18°. 
In order to obtain separated mode travel times, spectrograms were computed for various 
SUS explosive data. The theoretical mode travel times, for the ideal waveguide were then 
computed as a function of range R between source (SUS) and receiver (VLA), water depth 
H, water sound speed c,, and critical angle 0... A comparison between theoretical and 
measured travel times for all 11 SUS data sets produced an average value for the density 
ratio of 1.6. However, the changes in the theoretical curves for various density ratios was 
not large whereas the mode arrival structures were fairly ambiguous. Thus the uncertainty 
is very large. It is also possible that the theoretical mode travel times have significant errors 
due to the mismatch in sound speed profile. However, the estimated value of p, ~ 1.6 g/cm’ 
from the inversion method seems quite realistic. The results supplied by Miller and Potty 
suggest a more reasonable estimate for the bottom density is about 1.8 g/cm’. Figures 9a 
and 9b show the mode travel times for a phone from a SUS explosive data and those for an 


ideal waveguide. 
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c. INVERSION OF BOTTOM ATTENUATION 

A method was developed in Section III.C. to invert for bottom attenuation using 
measures of modal transmission loss as a function of range. In order for this method to 
work successfully, the decrease in amplitude due to bottom attenuation must be observable 
over the transmission loss due to cylindrical spreading. For example, as the range between 
the SUS drop location and the VLA increases from 15 km to 30 km, the change in 7L due 
to cylindrical spreading is approximately 3 dB. This equates to about 0.2 dB/km. In order 
for bottom attenuation to be distinguishable, Eq. (38) suggests that the quantity 20N ¢ 
should be of this same order of magnitude or larger. 

A simple analysis assuming bottom attenuation increase linearly with frequency 
shows that € increases with increasing frequency. Specific values can be obtained using 
Eqs. (34) - (44) with realistic values of bottom attenuation a. However, as frequency 
increases, the angle of propagation for a specific mode decreases thereby decreasing N, the 
number of bottom interactions per km. 

If one assumes a reasonable value of a ~ 0.05 dB/km/Hz, a 100 Hz signal suffers a 
0.5 dB/km loss due to attenuation for angles of propagation of about 12°. At this frequency 
in this waveguide, that corresponds roughly to Mode 2. Unfortunately, this mode is not 
clearly distinguishable at this frequency from the available SUS data. The lowest modes 
can be isolated at lower frequencies, but then the attenuation is negligible. At higher 
frequencies, the propagation angles for these low modes is near grazing, and again the 
influence of attenuation is small. 

The result of this analysis, unfortunately, shows that this data does not avail itself to 
inversion for bottom attenuation using the method developed previously. Other processing 
methods may be able to extract the necessary information in order to utilize this method, but 


will not be explored further here. 
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Figure 9a. Mode Travel Times for a Phone from a SUS Explosive Data. 
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Figure 9b. Mode Travel Times for a Phone from an Ideal Waveguide. 
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V. CONCLUSION 


This thesis incorporated several methods which may be applied directly to measured 


data to provide a relatively simple means of obtaining reasonable estimates of ocean bottom 


parameters from minimal information. The analysis of the data from the Mid-Atlantic 


Bight experiment recorded on the NW VLA produced the following results: 


l. 


Ambient noise was used to measure both the vertical directionality and spatial 
coherence of the noise field. By examining the noise data, we have been able to 
estimate 8. and hence the sound speed c, in the sediment. The bottom sound 
speed deduced was approximately 1563+10 m/s. 

The method of effective depth was employed to estimate the bottom density p,,. 
Although this method was derived for an isospeed water column environment, it 
is expected to produce reasonable estimates for most shallow water 
environments. By comparing mode travel times extracted from SUS data 
spectrograms and comparing these with values for an ideal waveguide with an 
effective depth, the bottom density was estimated to be p, ~ 1.6 g/cm’. The 
uncertainty in this value is very large, however, partly due to the difficulty in 
distinguishing mode arrivals clearly in the spectrograms. 

In order to distinguish bottom attenuation from simple cylindrical spreading, it 
was determined that values for mode transmission loss were required for high 
modes at high frequencies. Unfortunately, the SUS spectrogram data could not 
adequately provide this information, and no estimate of bottom attenuation was 
produced. It may be possible to use the method outlined if a modal 
decomposition of the data on the VLA could be performed. However, due to the 
relatively short aperture of the array (compared to the waveguide depth), this 
seems unlikely. Furthermore, the source level of different SUS detonations may 
not be repeatable enough to use this data for accurate estimations of bottom 


attenuation. 
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Comparing the inversion results for bottom sound speed and density with values 
obtained from genetic algorithm techniques provided by the University of Rhode Island, 
good agreement was found particularly for the sound speed. However, the analysis of the 
SUS data proved more difficult than originally anticipated. Much of this has to do with the 
uncertainty in absolute drop location and repeatable source signals. Another part of the 
problem was the reliance on low mode information which is difficult to characterize for 
SUS data measured on only eight hydrophones with total aperture less than half the width 
of the waveguide. It is likely that this analysis would be more successful with different 


more well-defined sources and a full water column array. 
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